arXiv:1508.03440vl [quant-ph] 14 Aug 2015 


Effects of electric field polarizations on pair production 


Z. L. Li/ D. Lu/ and B. S. 

^Key Laboratory of Beam Technology and Materials 
Modification of the Ministry of Education, 

College of Nuclear Science and Technology, 

Beijing Normal University, Beijing 100875, China 
^Beijing Radiation Center, Beijing 100875, China 
(Dated: August 17, 2015) 

Abstract 

The effects of electric field polarizations on pair production from a vacuum are investigated nu¬ 
merically by employing the real-time Dirac-Heisenberg-Wigner formalism. For few-cycle fields, it 
is found that the interference pattern in momentum spectra is absent and the circular distortion 
of momentum distribution becomes more apparent with the increase of polarization. For multi¬ 
cycle fields, it is found that the interference effects in momentum spectra are obvious. And as 
the polarization increases, the momentum distribution is split into two parts in the momentum 
Qy direction first and then two separated parts are connected into a ring. The effects of polar¬ 
izations on the number density of created particles exhibits two different characteristics. For a 
small laser frequency, the particle number density deceases with the polarization, while for a large 
laser frequency, the relation between them is sensitive to the field frequency nonlinearly. Some 
interpretations for the obtained results in this study, which is expected to be valuable to deepen 
the understanding of pair production in complex fields and also helpful for the study associated to 
strong-held ionization. 
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I. INTRODUCTION 


Electron-positron (EP) pair production from vacuum in the presence of a strong field is 
one of the well-known predictions of quantum electrodynamics (QED) [IH3] as a valuable 
probing the structure of QED vacuum. According to Schwinger’s formula, the pair produc¬ 
tion rate (oc exp(— vrEcr/E)) is exponentially suppressed for a weak held. So it seems that the 
experimental detection of vacuum pair creation is rather difficult because current electric held 
strength E is far less than the critical electric held strength E^r = vn? je 1.32 X 10^®V/cm 
(the units h = c = 1 are used). However, the high-intensity and ultrashort laser facilities 
under construction, such as the Extreme Light Infrastructure |1] and the x-ray free electron 
laser (XFEL) [5l E], may provide a subcritical electric held which can produce a consider¬ 
able number of EP pairs Ei. This stimulates the researcher renewed interests to explore 
signihcant phenomena of EP pair production from vacuum BM- 

In recent years, many theoretical works have been done about the pair production in 
more complex and realistic helds. These studies can not only be available to the future 
experiment but also reveal many novel signatures which are useful for understanding the 
complicated physics in pair production. Among existed studies the polarization of helds is 
a key parameter. Let us recall some important results. 

For a linearly polarized electric held, the momentum spectrum of created EP pairs in 
laser pulses with subcycle structures are investigated in Ref. [IS] with the help of quantum 
kinetic approach. It is found that the longitudinal momentum spectrum is extremely sen¬ 
sitive to the laser held parameters and its oscillations can be used as a probe of subcycle 
strucuturs in ultrashort laser pulses. Using the same electric held and calculation method, 
in Ref. [16], EP pair production in nonperturbative multiphoton regime are studied and 
the physical observable that carried the ehective mass signatures were identihed. For a cir¬ 
cularly polarized electric held that solely changes over time, Schwinger pair production is 
explored numerically with the real-time Dirac-Heisenberg-Wigner (DHW) formalism in Ref. 
la. and it is found that the characteristic momentum distribution in this held might be 
used to distinguish a QED cascade seeded by pair production from that seeded by isotropic 
vacuum impurity. Moreover, the pair creation rate for rotating electric helds is investigated 
semiclassically with the Wentzel-Kramers-Brillouin-type (WKB) approximation in Ref. [T8] . 
It is found that for a certain range of electric held parameters the pair production rate was 
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dominated by one spin of created pairs. 

Let us turn to the case, i.e., the problem of EP pair production in a general elliptically 
polarized electric held. In fact we have performed some primary studies a few years ago. 
For example, we studied analytically the exponential factor of pair production by using the 
worldline instanton method PI and a two-level transition technique PI. respectively. On 
the other hand, the elliptically polarized laser helds are also commonly used to explore new 
signatures of strong-held ionization (SFI) dynamics that can not be revealed by using a 
linearly polarized held [^Tl - E^ . Due to the similarity between strong-held ionization and 
vacuum pair production, the study of EP pair creation in an elliptically polarized held 
can not only uncover some new phenomena of pair production but also be a reference 
for strong-held ionization. Furthermore, as a elliptically polarized laser held is easier to 
achieve experimentally than a perfect circularly polarized one, the former is more suitable for 
clarifying pair production mechanism as well as observing EP pair creation experimentally. 

Motivated by the facts mentioned above that the importance of general elliptic polariza¬ 
tion exists in both of EP pair creation and SFI, therefore, in this paper, we focus our study 
on the ehects of electric held polarizations on EP pair production by numerically solving the 
real-time DHW formalism [25H27]. We calculate the momentum spectrum and the number 
density of created particles in diherent polarized electric helds with a same laser intensity 
and consider the ehects of held polarizations on them. We hnd that with increasing polariza¬ 
tions the circular distortion of momentum distribution becomes more obvious for few-cycle 
helds while the momentum spectra are split for many-cycle helds. It is also found that the 
particle number density deceases with the increase of polarization for a small laser frequency 
and the relation between them is more complicated for a large laser frequency. These results 
are helpful for the understanding of pair creation in complex helds. Evenly the present study 
may also be used to build a bridge between vacuum pair creation and SFI. 

This paper is organized as follows. In Sec. |^we give a brief recall of the DHW formalism 
which is used in our calculation. The numerical results are presented and analyzed in Sec. 
mil Section [13 is the conclusion and discussion. 
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II. THEORETICAL MODEL AND FORMULAS 


For a subcritical electric field E ~ O.li^cr; it is a good approximation to neglect the 
collision effect and the internal electric held since the EP pair yield and the back-reaction 
electric current are quite small. And because the spatial scales of the EP pair production 
are smaller than the spatial focusing scales of the laser pulse, the spatial effects are not 
signihcant. Therefore, we have the spatially homogeneous and time-dependent helds. For 
our studies, we focus on the EP pair production in a uniform and time-varying electric held 
of arbitrary polarization 
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Eq 


exp 



COs{ujt + 0 ) 

5 sin(a;t -|- 0) , 
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( 1 ) 


where Eq is the maximal held strength, r dehnes the pulse duration, oj is the laser frequency, 
0 is the carrier-envelope phase, and |5| < 1 represents the polarization (or the ellipticity). 
Note that a similar low-intensity electric held with the polarization up to ±0.93 has been 
achieved experimentally |2I]. For convenience, we set 0 = 0 and h > 0 unless otherwise 
specihed. 

Our following numerical results are based on the DHW formalism which has been used 
to study vacuum pair production in Refs. [25l |26] for diherent electric helds. We start with 
the equal-time density operator of two Dirac held operators in the Heisenberg picture. 


C«;3(x,y,t) = e--/-f/2A{x+Ay,q.ydA 


X 


x± -,t X- -,t 


( 2 ) 


with the center-of-mass coordinate x = (xi ±X 2)/2 and the relative coordinate y = xi — X 2 . 
Note that the factor before the commutator is a Wilson-line factor used to keep gauge 
invariance, and the integration path of the vector potential A is a straight line chosen to 
introduce a clearly dehned kinetic momentum p . Moreover, we have employed a Hartree 
approximation for the electromagnetic held and chosen the temporal gauge Aq = 0. The 
Wigner operator is dehned as the Fourier transformation of Eq. ([^ with respect to the 
relative coordinate y, and its vacuum expectation value gives the Wigner function 

>V(x,p,t) = -^ j dV*^'^(0|C(x,y,t)|0). (3) 
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Decomposing the Wigner function in terms of a complete basis set {1, 75, 7 ^, 7 ^ 75 , := 

I[ 7 ^, Y]}i we have 

>V(x,p,f) = ^(ls + i75p + 7 '"v^ + 7^75a + cr'"%j.), (4) 

with sixteen real Wigner components, scalar s(x, p,f), pseudoscalar p(x, p,f), vector 
v(x, p,t), axialvector a(x, p,t), and tensor t(x, p,t). Inserting the decomposition into the 
equation of motion for the Wigner funtion, one can obtain a partial differential equation 
(PDE) system for the sixteen Wigner components [25]. Furthermore, for the spatially ho¬ 
mogeneous and time-dependent electric helds mentioned above, by using the method of 
characteristics, or simply, replacing the kinetic momentum p by q — eA{t) with the well- 
dehned canonical momentum q, the PDE system for the sixteen Wigner components can be 
reduced to an ordinary differential equation system for the ten nontrivial Wigner components 
w(q, f) = (s,v,a,ti := 2 t*°ej)'''(q,f), 

w(q,t) = ?/(q,f)w(q,f), (5) 


where the dot denotes a total time derivative, 'H(q, f) is a 10 x 10 matrix. 
The one-particle distribution function is dehned as 

/(q.*) = ■ |w(q.i) - Wvaclq,i)l. 


( 6 ) 


where Wvac(q)^) (®vac ^vac 0 0) , ffivac 2?7l/D(p) |p_>q_eA(q) ^vac 

-2p/D(p)|p^q_eA(t), D(p)|p^q_eA(t) = + [q ” eA(f)]2)V2 jg the total energy of elec¬ 

trons, and ei = —1/2 Wvac is one of the basis of the ten-component vector w. Notice that 
the vacuum solution is Wvac(q 5 ^vac)- 

In order to precisely obtain the distribution function /, we employ the method used in 
m- Decomposing the Wigner components as w = 2(/ — 1 )(Ei -|-J^wg with an auxiliary nine- 

-p'''/m 0 


component vector Wg and a 10 x 9 matrix T = 
(|^, we have 


p^q—eA(t) 


, and applying Eq. 


/ = 1/2 (e[J^W9, 

Wg = 'HgWg 2(1 - /)^<Ei, 


( 7 ) 
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where ^ = (0 I 9 ) is a 9 x 10 matrix, and 


•Hg 


^ —ep ■ E''"/a;^(p) —2px —2m ^ 

-2px 0 0 

^ 2(m^ + p ■ p''')/m 0 ^ / 


p^q—eA(t) 


Thus, we can get the one-particle momentum distribution function /(q, t) by solving Eq. 
Q with the initial conditions /(q, —cxd) = W 9 (q, —cx)) = 0. Integrating the distribution 
function over full momenta at f —)■ -f-oo, we have the number density of created pairs 
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III. RESULTS 

Before analyzing our numerical results, we introduce the Keldysh adiabatic parameter 
7 = mu/eE [28] which is used to describe the dominant mechanism in pair production or 
strong-held ionization. For 7 1, the tunneling process dominate pair production and it 

is called nonperturbative Schwinger pair production. For 7 3> 1, the main contribution 
to pair creation comes from multiphoton mechnism and this process is called perturbative 
multiphoton pair creation. However, when 7 ~ 0{1), multiphoton pair production will have 
a nonperturbative nature and many interesting features of pair production will be presented 
[T6l [29l l30] . This intermediate regime is also called nonperturbative multiphoton regime. 

By solving Eq. ([^ numerically, we obtain the momentum spectra of created particles 
in the polarization plane Qy) for different polarizations with 7 = 0.35a/1 -|- 5^, see Fig. 

From Fig. [^a), one can see clearly the oscillatory structures of momentum spectrum. 
This result is similar to the one in P!. because for h = 0 our electric held ([^ will become 
the linearly polarized held used in [15]. When 5 7 ^ 0, however, the oscillations vanish and 
the peak position of momentum distribution is shifted along the direction of qy, see Fig. 
Kb)- With the increase of polarization 5, the momentum distribution of created EP pairs 
is distorted gradually. For h ~ 1, it forms a ring-like structure. For this few-cycle held 
(cr = car = 5), the absence of interference ehects in momentum distribution is because the 
polarization can easily change the distribution of tuning points (complex time tp satisfying 
= 0 ) cind as a result only one pair of turning points dominates the pair creation. 
Moreover, since the electric held component Ey{t) is an odd function of time, the momentum 
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distribution is not symmetric about Qy = 0 [3T], this can also be seen in Fig. 4 of Ref. [T5] . 
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FIG. 1: (color online) Momentum spectra of created EP pairs in the {qx,Qy) plane for different 
polarizations with to = 0.05m. (a) for 5 = 0, (b) for 5 = 0.5, (c) for 6 = 0.9, and (d) for 5 = 1. 
Other electric field parameters are chosen as Eq = {).l^/2Ecr and r = 100/m. 

Based on the analysis in Ref. na, it is known that for a circularly polarized field the cir¬ 
cular distortion of momentum distribution can be enhanced by increasing the pulse duration, 
because the produced particles have enough time to follow the rotation of this field. And 
when the circularly distorted momentum distribution closes into a ring, interference effects 
will occur. An intuitively plausible understanding about this is that the quantum wave func¬ 
tions of left-moving and right-moving particles carrying different phases will superpose each 
other and induce the quantum interference. In fact, this effect can not be fully understood 
just by using circularly polarized fields. To solve this problem, we increase the cycle number 
of electric field {a = 15) by increasing the pulse duration and plot the momentum spectra of 
created pairs for different polarizations in Fig. Indeed, as discussed above, we can see a 
ring-like structure and the interference effect in the momentum spectrum for 5 = 1, see Fig. 
[^c). Meanwhile, however, we also find that for a small value of 5 the particle distribution in 
momentum space is split into two segments in the qy direction and interference effects exist 
in each segment. This result indicate clearly that the interference effect is related to the 
cycle number of electric field. Actually, it is caused by the interference among different pair- 
creation amplitudes corresponding to different turning points in the complex t plane [T8] . 
For a many-cycle field, there will be many pairs of turning points having nearly the same 
distance to the real t axis and this can induce evident interference [HI]. Consequently, the 
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interference effect in momentum spectrum will become more obvious with the cycle number 
of electric held increasing. 

To see the effect of polarization on momentum spectra clearly, we show the momentum 
distribution of created pairs scanning over polarization with 7 = 0.71a/1 + in Fig. One 
can see that as the value of polarization increases the momentum distribution is split into 
two parts gradually in the qy direction, and then these two separated parts are connected 
to form a ring. More notably, when the value of polarization 5^3, the created particles 
with small momentum vanish. These hndings are very similar to the results of strong-held 
ionization of helium using an elliptically polarized laser pulses (cf. Fig. 1 in Ref. [22] )■ 
Here we give an intuitive understanding about the split of momentum distribution shown 
in Fig. 1^ According to Schwinger’s work |3], it is known that the pair production is 
suppressed exponentially for a weak held E due to the pair creation rate ~ exp^—nEcr/E). 
Therefore, for small values of polarization S, pair production manly occurs near the maximum 
of the electric held component Ex{t) and the electric held component Ey{t) plays the role 
in accelerating created particles. Because of the phase diherence between E^it) and Ey(t), 
the particles created at the peak of the hrst half cycle of E^it) are accelerated by Ey(t) 
along the y direction, while the particles created at the nadir of the second half cycle 
of Exit) are accelerated by Ey(t) along the —y direction. As a consequence, the created 
particles at diherent positions will have diherent momentum component qy, which results 
in two peaks of the hnal momentum distribution. This also explains why the distance of 
these two peaks increases with the value of polarization 6 increasing. For large values of 
polarization, the particles created by Ey{t) can also be split by the acceleration of Ex{t), 
thus the previous two separated parts are connected again. Moreover, for a few-cycle field, 
pair production is dominated by the held component Ex{t) at t = 0, so there is no split of 
momentum distribution, see Fig. For a many-cycle held, however, near t = 0 the electric 
held component Ex{t) has several cycles having almost the same amplitude as F^ 3 ;( 0 ), which 
induces the split of momentum distribution, see Figs. [^and|^ 

Generally, EP pairs created from vacuum via tunneling process are considered to be at 
rest and then are accelerated by strong laser pulses. However, on the basis of the similarity 
between pair creation and strong-held ionization (cf. Fig. [^here and Fig. 1 in Ref. [22]). 
the particles created from vacuum may have nonzero momentum. Of course, further studies 
are needed to prove this conjecture. Since many experiments about strong-held ionization 
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FIG. 2: (color online) Momentum spectra of created EP pairs in the {qx,Qy) plane for different 
polarizations with to = 0.05m and r = 300/m. From left to right 6 = 0.5, 0.9, and 1, respectively. 
The electric field strength Eq = 0.1\/2Ecr. 
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FIG. 3: (color online) Momentum spectra of created EP pairs in the {qx,Qy) plane for different 
polarizations with a; = 0.1m. The polarization value is shown in the headline of each panel. Other 
electric field parameters are the same as in Fig. 

are also carried out at 7 ~ 0{1), some explorations about nonperturbative multiphoton 
regime are useful. We know that in this intermediate regime both nonperturbative and 
perturbative process play an role in pair creation, so neither of them can fully explain the 
phenomena occurring in this regime. Thus a theory combining nonperturbative process with 
perturbative process will be hopeful to understand the experiment results |22lEl]. 

The momentum distributions of created EP pairs for different polarizations with 7 = 
3.54a/1 + are shown in Fig. For the large laser frequency u = 0.5m, one can clearly 
see some ring structures in the momentum spectra. These rings are formed by multiphoton 
absorbtion, i.e., absorbing n photons of frequency 00 to conquer the energy gap and produce 
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FIG. 4: (color online) Momentum spectra of created EP pairs in the {qx,Qy) plane for different 
polarizations with to = 0.5m. (a) for <5 = 0, (b) for 5 = 0.5, (c) for 6 = 0.9, and (d) for J = 1. 
Other electric held parameters are the same as in Fig. 

EP pairs. From the inside out, the rings correspond to 5-, 6-, 7-, 8-, and 9-photon absorp¬ 
tion, respectively. For the linearly polarized field (<5 = 0), we can calculate the ring radius of 
n-photon absorption qn = [{noo/2)‘^ — m^ — {m^Eo/ooEcr)'^ by employing the energy con¬ 
servation equation nu = 2 ml + with the effective mass m* = m[l -|- (mFo/a;Fcr)^/2]^/^ 
[TU] . For instance, for n = 5, the ring radius q^ ~ 0.7228m is in excellent agreement with 
the numerical result qnum ~ 0.7227m. Since we keep the laser intensity of the electric field 
([^ constant, the effective mass expression is approximately valid for arbitrary polarized 
helds according to the dehnition of effective mass m* = m[l -|- m?{—A^A ^)/with 
four-vector A^ [isiEg. So the ring radius can be roughly calculated by the expression for 
arbitrary polarized helds. For instance, the analytical result of the ring radius of 5-photon 
absorption in Fig. |^c) q^ ~ 0.7228m is approximately identical to the numerical result 
(Inum ~ 0.7267m. The small difference between these two results is caused by the effect of 
polarization which cannot be considered by the effective mass model. It is also found that 
with the value of polarization increasing the rings become more uniform and the created 
pairs with small momentum decreases due to the rotation of the electric held. 

From Figs. [T]j^ we can hnd that the maximum value of the momentum distribution func¬ 
tion f{+oc) decreases with increasing the value of polarization. So one might speculate that 
the number density of created EP pairs may also have the same change tendency. However, 
things are more complicated than we originally thought. To hgure out the relation between 
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FIG. 5: (color online) The number density of created particles n[qz = 0] as a function of the field 
frequency u for <5 = 0, 0.5, 0.9, and 1, respectively. Other electric field parameters are the same as 
in Fig. [T] 



FIG. 6: (color online) The number density of created particles n[qz = 0] as a function of the 
polarization 5 for different field frequencies. Other electric field parameters are the same as in Fig. 

m 

the number density and field polarization, we plot the number density scanning over the 
held frequency for different polarizations in Fig. It is found that there are many obvious 
oscillations on the curves and the ones corresponding to small held frequency vanish with 
the polarization increasing. These oscillations are caused by multiphoton pair production 
and have been studied for a linearly polarized held in Refs. Additionally, we also 

hnd that for a small held frequency the number density of created particles decreases greatly 
with the increase of the value of polarization. This is because for small laser frequencies non- 
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perturbative Schwinger pair production is more significant. As pair creation is suppressed 
exponentially for weak fields, the number density is dominated by the electric field compo¬ 
nent with a large amplitude. Therefore, the polarized field with a small value of S, which 
has a larger amplitude of field component than the one with a large value of 6, gives a high 
number density. For large field frequencies, we find that the relation between the number 
density and the polarization of electric field becomes more complicated because perturbative 
multiphoton process becomes more important. To see clearly, we plot the number density 
as a function of polarization for different laser frequency in Fig. It can be seen clearly 
that the number density of created pairs can reach its maximum at 5 = 0, 5 = ±1, or even 
0 < |5| < 1 for different field frequencies. Thus, for a large laser frequency the relation 
between the particle number density and the polarization of electric field is sensitive to the 
field frequency. 


IV. CONCLUSION AND DISCUSSION 

In summary, we have investigated the effects of electric field polarizations on vacuum 
pair production numerically by employing the real-time DHW formalism. First, we consider 
the effect of polarization on the momentum spectra of created particles and find that for a 
multi-cycle field the momentum distribution is split into two segments in the qy direction 
and then connected into a ring with increasing the polarization. This result is similar to 
the one in SFI. By studying the effect of polarization on the number density of created EP 
pairs, we find that for a small laser frequency the particle number density deceases with 
the increase of field polarization. However, for a large laser frequency, the number density 
can reach its maximum value at linear polarization, circular polarization, or even elliptical 
polarization, which seems to be counterintuitive. 

Due to the use of a general polarized field, our studies not only reveal several charac¬ 
teristics of pair production which can not be captured by a linearly polarized or circularly 
polarized field, such as the split of momentum spectrum, but also link up pair produc¬ 
tion with the similar phenomena in SFI that some experiments are ongoing. Furthermore, 
because of the limitation in experiment instruments, a perfect circularly polarized field is 
harder to be produced than a elliptically polarized one. So our studies about the elliptic 
polarization are also more closely associated to the experiment interests. 
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